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ABSTRACT 

Multichannel  parametric  filters  are  currently  being  studied  as 
a  means  of  reducing  the  dimension  of  STAP  algorithms  for  in¬ 
terference  rejection  in  airborne  pulsed-Doppler  radar  systems. 
These  filters  are  attractive  to  use  due  to  the  low  computational 
cost  associated  with  their  implementation  as  well  as  their  near 
optimal  performance  with  a  small  amount  of  training  data  for 
a  stationary  environment.  However,  these  filters  do  not  perform 
well  in  certain  types  of  non-stationary  environments.  This  pa¬ 
per  presents  two  modifications  to  the  Space-Time  AutoRegressive 
(STAR)  filter  that  we  previously  proposed.  The  first  modification 
is  based  on  the  Extended  Sample  Matrix  Inversion  (ESMI)  tech¬ 
nique  and  is  used  in  the  presence  of  range  varying  clutter  which 
arises  from  the  use  of  non-linear  antenna  arrays  or  bistatic  radar 
systems.  The  second  modification  to  the  STAR  filter  is  for  use  in 
the  presence  of  hot  clutter  and  is  a  three-dimensional  STAP  al¬ 
gorithm.  Using  a  realistic  simulated  data  set  for  circular  array 
STAP,  we  show  that  the  modifications  to  the  STAR  filter  improve 
the  performance  when  in  the  presence  of  the  non-stationary  in¬ 
terference. 

1.  INTRODUCTION 

The  use  of  space-time  adaptive  processing  (STAP)  for  airborne 
radar  interference  mitigation  is  usually  limited  by  the  lack  of  sta¬ 
tionary  secondary  data  used  for  training  the  filter.  This  problem 
is  made  worse  when  the  radar  platform  is  operating  under  circum¬ 
stances  that  lead  to  additional  non-stationary  components  to  the 
interference.  Such  circumstances  include  the  use  of  a  non-linear 
or  non- side-looking  array  which  leads  to  a  range  variation  of  the 
clutter  statistics  or  the  presence  of  an  airborne  jamming  source 
which  leads  to  hot  clutter  or  terrain  scattered  interference. 

Partially  adaptive  STAP  filters  alleviate  this  problem  to  a  de¬ 
gree  by  taking  advantage  of  the  low-rank  nature  of  the  clutter. 
The  partially  adaptive  STAP  filters  use  fewer  degrees  of  freedom 
and  therefor  need  fewer  training  samples  than  the  fully  adaptive 
STAP  filter.  One  such  partially  adaptive  STAP  filter  that  is  dis¬ 
cussed  in  this  paper  is  the  Space-Time  AutoRegressive  (STAR) 
filter  [1].  The  partially  adaptive  STAP  filters  offer  an  improve¬ 
ment  over  the  fully  adaptive  STAP  filter  but  are  still  derived  based 
upon  the  assumption  that  the  interference  is  stationary.  When  the 
non-stationary  component  of  the  interference  follows  a  specified 
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model,  this  model  may  be  taken  into  account  to  derive  a  filter  to 
cancel  the  non-stationary  interference.  A  few  partially  adaptive 
STAP  algorithms  have  been  derived  to  account  for  range- varying 
interference  [2]  and  hot  clutter  [3]. 

Parametric  filters  (such  as  the  STAR  filter)  have  been  shown  to 
achieve  near  optimal  performance  with  a  small  amount  of  training 
data  when  the  interference  is  stationary  [4].  However,  the  perfor¬ 
mance  when  the  interference  is  non-stationary  leaves  much  room 
for  improvement.  In  this  paper,  two  extensions  of  the  STAR  fil¬ 
ter  to  account  for  both  range- varying  interference  and  hot  clutter 
are  presented.  The  improvements  that  the  range- varying  Extended 
STAR  (ESTAR)  filter  offers  over  the  standard  STAR  filter  is  illus¬ 
trated  with  a  synthetic  data  set  generated  by  MIT  Lincoln  Labora¬ 
tory  that  simulates  the  output  of  a  20  element  antenna  array  whose 
elements  lie  along  a  circular  arc  of  120°  [2].  This  ESTAR  filter 
is  also  shown  to  have  better  performance  than  a  range- varying  ex¬ 
tended  post-Doppler  algorithm. 

The  three-dimensional  STAR  filter  used  to  mitigate  hot  clutter 
is  tested  using  the  same  data  set  as  above  augmented  with  syn¬ 
thetic  hot  clutter.  The  3D-STAR  filter  achieves  a  significant  im¬ 
provement  in  signal-to-interference  plus  noise  ratio  (SINR)  over 
the  standard  STAR  approach.  In  comparing  the  3D-STAR  filter  to 
a  three-dimensional  optimized  pre-Doppler  algorithm,  it  is  shown 
that  the  performance  of  the  two  filters  are  nearly  the  same  but  that 
the  3D-STAR  filter  has  a  narrower  clutter  notch.  This  narrow  clut¬ 
ter  notch  allows  for  improved  detection  of  slowly  moving  targets. 

In  the  next  section,  we  briefly  present  the  standard  data  model 
used  for  STAP  problems  and  introduce  the  notation  used  through¬ 
out  the  paper.  The  STAR  filtering  technique  is  described  in  Section 
3  as  a  background  for  the  extensions  presented  herein.  Section  4 
presents  the  range- varying  extended  STAR  filter  that  is  used  when 
the  clutter  statistics  are  range- varying.  Section  5  derives  a  3D- 
STAR  filter  used  for  the  mitigation  of  hot  clutter  and  Section  6 
shows  the  results  of  several  numerical  simulations  of  the  filters. 

2.  DATA  MODEL 

A  target  present  in  a  particular  range  bin  during  some  coherent  pro¬ 
cessing  interval  (CPI)  may  be  modeled  as  producing  the  following 
baseband  vector  signal  (after  pulse  compression  and  demodula¬ 
tion)  [5]: 

xrit)  =  6a(0)e^“'‘  +  nriri  e  cr ,  t  =  (1) 

where  i  is  the  range  bin  in  which  the  target  is  located,  h  is  the 
complex  amplitude  of  the  signal,  oj  is  the  Doppler  shift  due  to  the 
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relative  motion  between  the  array  platform  and  the  target,  a(^)  is 
the  response  of  the  array  to  a  unit  amplitude  plane  wave  arriv¬ 
ing  from  direetion  9  (azimuth  and  elevation  angles),  and  ni{t) 
contains  contributions  from  clutter,  jamming,  and  thermal  noise. 
In  (1),  we  are  assuming  an  array  of  m  elements  and  a  total  of  N 
transmitted  pulses  covering  R  range  bins. 

If  we  stack  the  N  array  outputs  into  a  single  mN  x  1  space- 
time  snapshot,  we  may  re-write  (1)  as 


Xi 


x(l) 

x(7V) 


=  bs{d,L0)  +  T] 


(2) 


where 


s{d,u>)  —  v(w)®a(0) 

v(w)  =  [1  •••  r 

ri  =  [n(lf  ...  n{Nff 

and  0  represents  the  Kronecker  product.  The  vector  rji  contains 
the  stacked  vector  samples  of  the  clutter  and  interference  for  range 
bin  i,  and  has  an  unknown  covariance  matrix  denoted  by 


In  cases  where  the  clutter  is  stationary,  we  assume  that  equations  (4) 
and  (5)  also  hold  for  the  secondary  data  as  well: 

H*N,  =  0  (7) 

=  0 ,  (8) 

for  A;  =  1,  •  •  • ,  A^5,  where  Ns  is  the  number  of  secondary  data 
snapshots  used  to  train  the  filter. 

The  matrix  is  mNxm'{N—L-\-l).  If  (3)  holds  andm'  and 
L  are  chosen  so  that  m'{N  —  L  +  1)  is  large  enough,  the  columns 
of  Ti  form  a  basis  for  the  space  orthogonal  to  the  clutter  and  in¬ 
terference  subspace.  Although  this  relationship  does  not  hold  in 
practice  due  to  the  presence  of  thermal  noise,  a  least  squares  so¬ 
lution  is  applied  to  approximate  the  subspace.  This  suggests  the 
following  space-time  filter  (similar  to  the  matched  subspace  detec¬ 
tors  in  [7])  be  used  for  interference  rejection: 

WA«(^,w)  =  ,  (9) 

where  P-h  is  the  projection  onto  the  columns  of  R: 

.  (10) 


=  R. . 

The  clutter  is  neither  temporally  nor  spatially  white;  in  fact,  the 
rank  of  R  is  typically  much  less  than  mN.  The  rank  (p)  of  R  is 
important  because  it  determines  how  many  secondary  data  sam¬ 
ples  are  required  to  accurately  estimate  R.  According  to  [6],  the 
number  of  required  samples  is  on  the  order  of  2/9  to  5/9.  The  fully 
adaptive  approach  to  whitening  this  type  of  data  is  to  multiply  the 
data  by  the  inverse  square-root  of  an  estimate  of  the  matrix  R.  Be¬ 
cause  the  size  of  this  matrix  can  become  quite  large,  its  low  rank 
nature  is  exploited  to  derive  reduced-dimension  whitening  algo¬ 
rithms.  The  next  section  summarizes  the  work  in  [1]  as  a  back¬ 
ground  for  extending  the  STAR  filter. 


3.  SPACE-TIME  AUTOREGRESSIVE  FILTERING 


Following  the  derivation  in  [1],  the  STAR  approach  assumes  that  a 
set  of  L  matrices  Ho,  Hi,  •  •  • ,  Hl-i  of  dimension  m'  x  m  exist 
that  satisfy 


L-l 

H^n(t  +  i)  =  0  ,  t  =  1,  -  ■  ■  ,N  -  L  1,  (3) 


for  the  interference  and  clutter  in  the  primary  range  bin.  We  may 
also  write  (3)  in  the  following  two  different  ways: 


[Ho  •••  Hl-i] 


1(1) 


i(L) 


n{N  -L  +  1) 
n{N) 


=  0  (4) 


N 


or 

n*ri  =  0, 

where 


n* 


Ho  Hl-1 

Ho  Hl-1 


Ho 


(5) 


(6) 


We  refer  to  the  implementation  of  STAR  with  the  weight  vector 
of  (9)  as  Space-Time  AutoRegressive  (STAR)  filtering.  The  STAR 
filter  weights  are  “adaptive”  in  the  sense  that  R  must  be  estimated 
from  the  secondary  data  prior  to  computation  of  war. 


4.  RANGE-VARYING  EXTENDED  STAR  FILTER 


The  STAR  filter  of  the  previous  section  is  not  designed  to  han¬ 
dle  non- stationary  interference  of  any  kind.  This  section  derives  a 
STAR-based  filter  that  assumes  the  clutter  statistics  vary  linearly 
with  range.  This  assumption  is  reasonable  if  the  training  region  is 
kept  short.  The  idea  of  using  time- varying  weights  in  a  STAR  algo¬ 
rithm  was  introduced  in  [8]  as  an  extended  sample  matrix  inversion 
algorithm  and  this  idea  was  used  for  range- varying  STAR  weights 
in  [2] .  This  technique  increases  the  dimension  of  the  problem  by 
a  factor  of  two  but  does  improve  the  performance  when  there  is  a 
rapidly  changing  clutter  locus. 

The  idea  behind  range- varying  weights  is  that  the  weight  vec¬ 
tor  is  a  function  of  range  (r)  to  account  for  the  non- stationary  clut¬ 
ter  locus.  Expanding  the  weight  vector  into  a  power  series  yields 

t‘^Wq 

w(r)  =  AVo  +  rwo  H - ^  - • 

The  assumption  is  made  that  the  clutter  locus  is  changing  slowly 
enough  that  for  a  given  collection  of  ranges  the  weight  vector  is 
linear  in  r.  Ignoring  the  higher  order  terms  in  the  Taylor  series, 
the  weight  vector  as  a  function  of  the  range  bin  becomes 

w/c  =  Wo  +  akAwo,  (12) 

where  a  is  a  normalization  constant.  Defining 


Wo 

A’w 


(13) 


the  output  of  the  filter  may  be  written  as 


Xk 

akxk 


Hl-1 


z  =  w* 


(14) 


where  x/c  is  the  extended  data  vector. 

Using  this  same  idea  for  the  STAR  filter  {i.e.,  assuming  that 
the  STAR  filter  coefficients  that  null  the  clutter  vary  linearly  with 
range)  we  can  rewrite  (3)  using  an  extended  data  vector: 


E[ 


Hi 


AHi  ] 


iik{t  -  i) 
akiik{t  —  i) 


=  0,  t  =  - 


(15) 

Letting  AH  and  be  defined  similar  to  H  in  (4)  and  l-i  in  (6) 
we  may  rewrite  (7)  and  (8)  as 


[  H*  AH*  ] 

[  W  An*  ] 


N, 

aklSk 

rjk 

akrjk 


0  (16) 

0.  (17) 


The  filter  parameters  H  and  AH  may  then  be  estimated  using  the 
left  null  space  of  the  matrix 


N_g  •••  Ng 

-aQN_g  •  •  •  aQNg  J  ’ 


(18) 


where  Q  —  Following  what  was  done  in  [2],  the  constant  a 
is  chosen  as 

^ (TV,  +2) (TV,  ^ 

to  yield  a  “fiat”  noise  subspace. 

To  define  what  the  weight  vector  is,  let 


so  that 


^e{0,  (jj) 


s(^,  (jj) 

0 


(20) 


(21) 


Filtering  the  extended  data  vector  with  (21)  is  referred  to  as  the 
Extended  STAR  (ESTAR)  filter.  When  estimating  a  range  varying 
weight  vector  using  data  that  also  varies  with  range,  a  higher  num¬ 
ber  of  training  vectors  may  be  used  before  performance  starts  to 
degrade. 


5.  STAR  FILTERING  FOR  HOT  CLUTTER 

When  the  radar  platform  is  operating  in  an  environment  where 
there  is  an  airborne  jamming  source  present,  two  main  considera¬ 
tions  must  be  made.  First,  the  hot  clutter  covariance  changes  from 
pulse  to  pulse  and  second,  the  hot  clutter  has  non-zero  correlations 
across  range  bins  [3].  This  section  derives  a  STAR  based  filter 
that  is  effective  in  canceling  hot  clutter.  The  baseline  STAR  filter 
is  first  modified  to  handle  any  type  of  interference  that  changes 
from  pulse  to  pulse  (as  with  intrinsic  clutter  motion)  and  then  an 
additional  dimension  is  added  to  the  vector  autoregressive  filter  to 
account  for  the  correlations  across  range  bins. 

The  model  for  the  clutter  in  (3)  is  no  longer  valid  since  the  spa¬ 
tial  covariance  changes  from  pulse  to  pulse.  If  the  standard  STAR 
model  is  used  in  a  non- stationary  environment  like  hot  clutter,  it 
tries  to  account  for  the  time  variations  in  the  data  by  increasing 
the  number  of  filter  taps  required  to  achieve  clutter  cancelation. 


A  better  model  for  this  is  to  let  the  coefficients  of  the  space-time 
prediction  error  filter  change  with  time: 


r  Ho(l)  Hl-i(1) 


07*  _ 

ftTV  — 


Ho(n)  •••  HL-i(n) 


(22) 

where  n  —  N  —  L  -\-  1  and  where  each  block  row  is  a  set  of 
new  coefficients  based  on  dropping  the  data  from  the  oldest  pulse 
and  adding  the  data  from  the  most  recent  pulse.  For  this  time- 
varying  STAR  filter,  a  greater  number  of  filter  parameters  must  be 
estimated  {n  times  the  degrees  of  freedom  required  for  the  standard 
STAR  algorithm)  and  therefore,  more  sample  support  is  required 
to  train  the  filter. 


To  complete  the  derivation  of  the  3D-STAR  filter,  a  few  defi¬ 
nitions  need  to  be  made.  To  clarify  the  notation,  sampling  across 
pulses  is  called  slow-time  sampling  and  sampling  across  range 
bins  is  called  fast- time  sampling.  Let  P  be  the  number  of  fast¬ 
time  samples  over  which  the  hot  clutter  is  correlated. 

In  order  to  utilize  the  fast-time  correlation  of  the  data,  an  extra 
dimension  is  added  to  the  STAR  filter.  We  assume  for  a  moment 
that  the  interference  is  stationary  across  the  pulses  (slow-time). 
This  filter  will  model  the  fast-time  and  slow-time  correlations  with 
a  two-dimensional  VAR  filter.  For  a  set  of  L  J  matrices  of  size 
M'  X  M,  assume  that  the  clutter  obeys  the  model 


J-l L-l 

jnk+j{t-\-i)  =  0, 


j—O  i—O 


t  =  !,•••, W-L  +  1 
=  J  +  1,  (23) 


where  A;  =  0  is  the  range  bin  of  interest  and  n/e(t)  is  the  spatial 
snapshot  for  the  pulse  and  the  k^^  range  bin.  This  may  also  be 
expressed  as 


j-i 


^n*Gk+j 

j=0 


=  0  = 


P- J  +  1, 


(24) 


where  Pj  is  the  matrix  defined  in  (6)  with  a  subscript  j  to  indicate 
which  fast-time  sample  it  is  associated  with.  From  this  point  we 
again  take  into  account  the  slow-time  variations  caused  by  the  hot 
clutter  by  replacing  Pj  with  the  slow-time  varying  filter  Ptvj  • 
Rewriting  this  sum  with  the  time- varying  filter  we  get 


=0,  (25) 


where 


moik) 


r)k 

fjk+p-i 


ir 


‘U*  ‘U* 

rLTV,0  ■■■  I^TV,J-1 

^TV,0  ■■■  J 


Assuming  that  there  is  target  energy  in  the  A;  =  0  range  bin,  then 
there  will  also  be  target  energy  in  the  vectors  773^(0),  773d(— 1), 
•  •  •,  f]^d{—P  +  1)  which  may  not  be  used  for  training  the  filter.  In 


order  to  define  the  algorithm  to  find  the  filter  coefficients  let 


H(t)*  =  [  Ho,o(t)  ]  (26) 

nfc(t) 

Skit)  =  :  (27) 

nfc(^  +  L  —  1) 

Skit)  ■■■  Sk+p-jit) 

Gkit)  =  :  :  (28) 

.  Sk+J-lit)  •••  Sk+P-lit)  _ 

Git)  =  [  ai(t)  •••  aiv,(t)].  (29) 


The  filter  coefficients  can  then  be  found  by  the  following  least 
squares  criterion: 

H(t)  =  argmin  ||H(t)*G(t)||  t  =  — L-\-l  (30) 

H(t)  II  IIf 

subject  to  the  constraint  that  =  I.  From  this  point 

the  m'  left  singular  vectors  corresponding  to  the  smallest  singular 
values  of  each  G{t)  matrix  will  be  used  to  compute  the  —  L  + 1 
sets  of  filter  coefficients  which  define  HL  With  a  defined  sub  space, 
a  weight  vector  for  mitigation  of  hot  clutter  is 

W3£)(^,6i;)  =  PHS3D(^,a;),  (31) 

where 

■  1  ■ 

0 

S3d{0,(jj)=  .  (32) 

_  0  _ 

This  3D-STAR  filter  will  require  more  training  data  than  the  STAR 
filter  (on  the  order  of  W  —  L  +  1  times  more)  due  to  the  non¬ 
stationary  prediction  error  filter  that  is  used  in  the  implementation. 
This  additional  sample  support  requirement  is  less  of  an  issue  than 
with  other  3D  implementations  because  the  STAR  approach  typ¬ 
ically  requires  much  less  secondary  data  for  good  performance. 
The  3D-STAR  filter  also  assumes  that  the  data  is  stationary  for  P 
fast- time  samples. 

6.  NUMERICAL  RESULTS 

The  algorithms  presented  herein  are  tested  using  a  data  set  created 
by  MIT  Lincoln  Laboratory  that  simulates  the  output  of  a  20  ele¬ 
ment  array.  These  elements  lie  along  a  circular  arc  of  120°  with 
radius  2.96m  and  are  assumed  to  have  a  cosine-shaped  response 
with  a  -30  dB  backlobe  for  both  azimuth  and  elevation  dimensions. 
The  airborne  platform  is  moving  with  a  velocity  of  100  m/s  above 
a  4/3  eath  model  at  an  altitude  of  9000m.  The  operating  frequency 
of  the  radar  is  taken  to  be  435  MHz,  the  radar  bandwidth  and  sam¬ 
pling  frequency  are  3.75  MHz,  the  pulse-repetition  frequency  is 
300  Hz,  and  N  =  IS  pulses  are  assumed  to  be  transmitted  dur¬ 
ing  one  CPI.  Data  are  generated  for  9325  range  gates  between  20 
and  400  km  with  a  clutter-to-white-noise  power  ratio  of  40  dB  at  a 
range  of  100km. 

Hot  clutter  is  included  in  the  data  by  adding  a  term  of  the  form 

C^l) 

ik  =  bj  : 

.  <^kiN)  . 


where  bj  is  the  amplitude  of  the  jammer, 

i 

Ckit)  =  a{0j)zk -\-'^h{t)zk-i 
i—l 

is  the  contribution  of  the  hot  clutter  for  a  single  pulse  at  range  k, 
i  is  the  longest  multipath  delay,  Oj  is  the  direction  of  arrival  of  the 
jammer  signal,  Zk  is  the  jammer  waveform  (white  in  both  slow  and 
fast- time),  and  bt  is  a  random  vector  that  approximates  the  sum  of 
the  spatial  steering  vectors  for  each  of  the  multipath  signals.  When 
present,  the  jammer- to-clutter  power  ratio  is  assumed  to  be  10  dB. 
When  secondary  data  are  used  to  estimate  the  clutter  covariance  or 
STAR  filter  parameters,  equal  amounts  of  data  from  range  bins  on 
either  side  of  the  target  range  bin  are  used. 

The  true  clutter  covariance  matrix  used  to  generate  the  data 
is  known  for  20  of  the  9325  range  bins,  and  thus  the  maximum 
achievable  SINR  can  be  calculated  at  these  ranges.  To  illustrate 
the  performance  of  the  algorithms  we  use  either  the  SINR  loss 
as  a  function  of  Doppler  for  an  azimuth  of  0°  or  the  “average” 
SINR  loss  as  compared  with  the  optimal  (known  covariance)  so¬ 
lution.  This  average  SINR  loss  is  defined  as  the  area  between  the 
algorithm’s  SINR  curve  and  that  achievable  assuming  R  is  known. 
This  is  depicted  in  Figure  1 .  The  ESTAR  filter  will  be  compared 
to  the  range- varying  extended  post-Doppler  PRI  staggered  STAP 
algorithm  [2]  and  the  3D-STAR  algorithms  will  be  compared  with 
the  optimized  3D  pre-Doppler  STAP  algorithm  [3].  For  the  STAR 
based  filters,  M'  —  20  is  used  for  all  the  examples  and  for  the 
partially  adaptive  STAP  algorithm,  three  pulses  at  a  time  are  pro¬ 
cessed  and  a  diagonal  loading  of  about  five  times  the  noise  level  is 
used  for  sample  matrix  inversion. 

A  performance  evaluation  of  the  ESTAR  filter  at  a  range  of 
20km  is  shown  in  Eigures  2  and  3.  Eigure  2  compares  the  perfor¬ 
mance  of  the  ESTAR  filter  and  the  basic  STAR  filter  as  a  function 
of  L  for  Ns  —  50  (2km  training  window).  This  figure  shows 
that  the  ESTAR  filter  does  perform  better  than  the  STAR  filter  at 
close  ranges.  We  also  see  that  the  ESTAR  filter  requires  fewer  fil¬ 
ter  taps  than  the  STAR  filter  thus  offsetting  some  of  the  additional 
computational  cost  associated  with  the  extended  implementation. 
Eigure  3  compares  the  performance  of  the  STAR  filters  with  the 
range-varying  extended  PRI  staggered  and  fully  adaptive  STAP 
algorithms  as  a  function  of  training  length.  Note  that  the  perfor¬ 
mance  of  the  STAR  algorithm  degrades  quickly  as  more  training 
data  is  used.  The  extended  PRI  STAP  and  ESTAR  filters  both  have 
nearly  fiat  performance  as  Ns  is  increased  due  to  the  range- varying 
weights.  The  ESTAR  filter  also  has  much  better  performance  than 
the  extended  PRI  STAP  algorithm  because  it  requires  much  less 
training  data  to  converge  to  its  best  performance. 

Another  aspect  of  performance  is  the  computational  load  re¬ 
quired  to  implement  the  algorithms.  Eor  the  STAR  algorithms  the 
implementation  is  broken  up  into  two  steps.  The  first  step  involves 
taking  the  SVD  of  the  2 ML  x  {N  —  L  l)Ns  data  matrix  M 
and  the  second  is  forming  the  projection  operator.  The  bulk  of  the 
computation  involved  in  this  second  step  is  finding  the  inverse  of 
iCl-l  which  is  usually  a  sparse  banded  matrix.  Taking  this  into 
account  the  computational  load  for  the  ESTAR  algorithm  is 

0(4(ML)^(7V  -  L  +  \)Ns)  +  0((ML)^(iV  -L  +  \)M'). 

For  the  parameters  of  the  circular  array  data  with  L  =  4  and  M'  = 
20,  the  computational  cost  is 

ESTAR  =  0(3.84  x  10^ Ns)  +  0(1.92  x  10®). 


Comparing  this  with  the  cost  of  the  STAR  filter  (at  L  =  5): 

STAR  =  0((ML)\N  -  L  +  l)Ns) 

+OiiML)^{N  -L  +  1)M') 

=  0(1.4  X  10® A,)  +  0(2.8  X  10®) 

the  ESTAR  algorithm  has  only  a  small  increase  in  computational 
load.  The  extended  PRI  STAP  algorithm  has  a  computational  cost 
of 

EPRISTAP  =  0{4{MK)^{N -K +  l)Ns) 

+0(4(MKf(N  -K  +  l)p) 

=  0(2.3  X  10® A,)  +  0(2.0  X  lo’’) 

where  K  =  3  pulses  that  are  processed  at  a  time  and  p  =  90  is  the 
approximate  rank  of  each  sub-CPI.  From  this  we  see  that  if  Ng  is 
not  too  big  (Ns  <  100),  then  the  ESTAR  algorithm  requires  much 
fewer  computations  than  the  PRI- staggered  STAP  algorithm. 

Figures  4-6  illustrate  the  performance  of  the  3D-STAR  filter 
when  there  is  hot  clutter  present  and  when  the  direct  path  jamming 
signal  is  in  the  mainbeam  of  the  radar  system.  Figure  4  compares 
the  performance  of  the  3D-STAR  filter  to  the  basic  STAR  filter 
as  a  function  of  L.  The  3D-STAR  filter  outperforms  the  STAR 
filter  with  a  small  number  of  filter  taps  by  utilizing  the  slow-time- 
varying  taps  as  well  as  the  additional  fast  time  tap.  Figure  5  com¬ 
pares  the  STAR  filters  to  the  3D  optimized  pre-Doppler  and  fully 
adaptive  STAP  algorithms  as  a  function  of  training  data.  In  this 
case  the  pre-Doppler  and  3D-STAR  algorithms  have  a  very  similar 
performance  with  the  pre-Doppler  algorithm  slightly  outperform¬ 
ing  the  3D-STAR  filter.  However,  Figure  6,  which  shows  the  SINR 
at  Ns  =  80  or  3.2  km,  illustrates  that  the  3D-STAR  filter  has  a 
narrower  clutter  notch  which  results  in  a  lower  detectable  veloc¬ 
ity.  If  the  small  loss  in  performance  away  from  the  clutter  notch 
is  tolerable,  the  3D-STAR  filter  is  more  desirable  due  to  its  greater 
percentage  of  usable  Doppler  space. 

The  computational  cost  of  the  STAR  (L  =  7),  3D-STAR  (L  = 
2,  J  =  2),  and  3D-pre  Doppler  (K  —  3  pulses)  algorithms  for  the 
system  parameters  described  above  are  as  follows: 

STAR  =  0(2.35  x  10® A,)  +  0(4.7  x  10®) 

3D -STAR  =  0{{MLjf{N-L4-\){P-J  +  l)Ns) 

+0{{MLJ)\N  -L  +  1)(P  -  J  +  1)M') 
=  0(2.18  X  10® A,)  +  0(4.35  x  10®) 

pre-Dopp  =  0{{MKPf{N-K  +  l)Ns) 
+0{{MKPf{N  -K  +  l)p) 

=  0(5.18  X  10® A,)  +  0(7.0  X  10^) 

where  p  =  135  is  the  approximate  rank  of  the  sub-CPI  covariance 
matrix.  Again  we  see  that  the  STAR  and  3D-STAR  algorithms 
have  nearly  the  same  computational  cost  when  the  filter  orders 
are  chosen  close  to  the  best  value.  It  is  also  seen  that  the  pre- 
Doppler  algorithm  requires  a  large  number  of  computations  when 
compared  with  the  3D-STAR  algorithm. 

7.  CONCLUSIONS 


The  first  modified  filter  (ESTAR)  is  used  when  the  clutter  statistics 
are  varying  with  range  as  is  the  case  for  non-linear  antenna  arrays 
or  bistatic  radar  systems.  The  second  modification  (3D-STAR)  is 
used  in  the  presence  of  hot  clutter  which  arises  when  an  airborne 
jamming  source  is  present.  These  two  modifications  provide  an 
increase  in  performance  over  the  standard  STAR  filter  when  used 
in  non- stationary  environments  without  a  major  increase  in  com¬ 
putational  burden.  We  have  shown  in  numerical  experiments  and 
computational  analysis  that  the  ESTAR  filter  is  superior  to  the  ex¬ 
tended  PRI-staggered  post-Doppler  STAP  algorithm  when  there  is 
a  rapidly  changing  clutter  locus.  We  have  also  shown  that  the  3D- 
STAR  filter  has  a  little  more  usable  Doppler  space  than  the  3D  op¬ 
timized  pre-Doppler  STAP  algorithm  and  the  3D-STAR  algorithm 
achieves  this  performance  with  much  less  computation. 

8.  REFERENCES 

[1]  A.  L.  Swindlehurst  and  P.  Parker,  “Parametric  clutter  rejection 
for  space-time  adaptive  processing,”  in  Proc.  of  the  ASAP- 
2000  Workshop,  (MIT  Lincoln  Laboratory,  Lexington,  MA), 
2000. 

[2]  M.  Zatman,  “Circular  array  STAP,”  IEEE  Trans,  on  Aerospace 
and  Electronic  Systems,  vol.  36,  pp.  510-517,  April  2000. 

[3]  D.  J.  Rabideau,  “Clutter  and  jammer  multipath  cancellation 
in  airborne  adaptive  radar,”  IEEE  Trans,  on  Aerospace  and 
Electronic  Systems,  vol.  36,  pp.  565-583,  April  2000. 

[4]  J.  Roman,  D.  Davis,  and  J.  Michels,  “Multichannel  parametric 
models  for  airborne  phased  array  clutter,”  in  Proc.  IEEE  Radar 
Conf,  vol.  1,  (Syracuse,  NY),  pp.  72-77,  1997. 

[5]  J.  Ward,  “Space-time  adaptive  processing  for  airborne  radar,” 
Tech.  Rep.  TR-1015,  MIT  Lincoln  Labs,  Dec  1994. 

[6]  I.  S.  Reed,  J.  D.  Mallett,  and  L.  E.  Brennan,  “Rapid  conver¬ 
gence  rate  in  adaprive  arrays,”  IEEE  Trans,  on  Aerospace  and 
Electronic  Systems,  vol.  AES- 10,  pp.  853-863,  Nov  1974. 

[7]  L.  Scharf  and  B.  Friedlander,  “Matched  subspace  detectors,” 
IEEE  Trans,  on  Sig.  Proc.,  vol.  42,  pp.  2146-2157,  August 
1994. 

[8]  S.  D.  Hayward,  “Adaptive  beamforming  for  rapidly  moving 
arrays,”  in  Int.  Conf  on  Radar,  (Beijing),  pp.  480-483,  CIE, 
Oct.  1996. 


This  paper  has  presented  modifications  to  the  space-time  autore¬ 
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Figure  3:  Comparison  with  PRI-Staggered  and  fully  adaptive  Figure  6:  SINR  of  3D-STAR  and  pre-Doppler  algorithms  with 
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